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ABSTRACT 

Many of the characteristic properties of the millisecond pulsars found in globular clusters are 
markedly different from those in the Galactic disc. We find that one such physical parameter is 
the surface magnetic field strength. Even though the average spin-periods do not differ much 
the average surface magnetic field is 2-5 times larger in the globular cluster pulsars. This 
effect could be apparent, arising due to one or more of several biases. Alternatively, if future 
observations confirm this effect to be real, then this could be interpreted as a preferential 
recycling of pulsars in tight binaries where the mass transfer takes place at high accretion 
rates. 

Key words: magnetic field - stars; neutron - pulsars: general - globular clusters: general. 



1 INTRODUCTION 



A radio pulsar is a strongly magnetized rotating neutron star. Dis- 
covered serendipitously (Hewish et al. 1968), the pulsars are char- 
acterized by their short spin-periods (P) and very large inferred 
surface magnetic fields (B). However, the detection of the 1. 5 ms 
pulsar B 1937+21 (Backer et al. 1982) heralded a new genre, that 
of the radio millisecond pulsars (MSP). The ranges of the spin- 
periods and the surface magnetic fields of these MSPs place them 
in a nearly disjoint region of the B — P plane from the normal radio 
pulsars. Close to two thousand radio pulsars have been detected to 
date with P ~ 10"^ - 10 s and B ~ 10* - 10^^ G. Amongst 
these, the MSPs are typically characterized by P< 30 ms, and typ- 



ical surface field strengths of ' 
spin-down ages of ^ 10^ yr. 
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The connection between these two seemingly disjoint populations 
of pulsars is realized through the binary association of the neutron 
stars. An MSP is understood to descend from an ordinary, long- 
period pulsar that have been spun-up and recycled back as an MSP 
in a low mass X-ray binary (LMXB) by mass accretion (Radhakr- 
ishnan & Srinivasan 1982). In the particular case of an isolated 
MSP, like PSR B1937+21, the donor is later destroyed probably by 
the wind of the newly recycled pulsar itself (Alpar et al. 1982). In 
general the MSPs are considered to be the end products of LMXBs 
and intermediate-mass X-ray binaries where the primary neutron 
star is assumed to have formed by the core-collapse of a mas- 
sive (Af > 8M0) star (Bisnovatyi-Kogan & Romberg 1974; Bhat- 
tacharya & van den Heuvel 1991). 

This scenario for MSP formation has been backed by several 
observational indications over the years. The strongest support 
for the connection between MSPs and X-ray binaries come from 
the discovery of coherent millisecond X-ray pulsations in SAX 
J1808.4 (Chakrabarty & Morgan 1998; Wijnands & van der Klis 



1998) and subsequently in many other systems (Wijnands et al. 
2005). These accreting millisecond X-ray pulsars (AMXP) are un- 
derstood to be the immediate precursors of the radio MSPs and are 
expected to turn into those when active mass accretion stops. The 
identification of PSR J1023-I-0038 as an MSP which has only re- 
cently lost its accretion disc (Archibald et al. 2009) provides almost 
a direct proof of this theory of MSP formation. 

However, it seems that the large number of pulsars detected in glob- 
ular clusters in the recent years have somewhat different character- 
istics than the pulsars observed in the Galactic disc. To begin with, 
the fraction of binary and millisecond pulsars in the clusters are 
~ 40% and ~ 92% respectively, whereas in the disc these fractions 
are only ~ 5% and ~ 4%. The size of the MSP population is much 
larger in the clusters than in the disc. In fact, the globular clusters 
(about 150 of them) known to orbit the Milky Way contain roughly 
three orders of magnitude more observed MSPs per unit mass than 
the Galactic plane, which contains 73 known MSPs (Camilo & Ra- 
sio 2005). The conditions prevailing in a globular cluster is rather 
different from that in the Galactic disc primarily due to the ex- 
tremely high stellar densities in the clusters. One of the obvious 
effect of this high density is a dramatic increase in the number of 
binaries, as well as in the rate of close stellar encounters allowing 
for many different channels for binary formation. Effectively, these 
systems provide different ways of pulsar recycling, than seen in the 
disc where the MSPs are mainly formed in primordial low-mass 
binaries. 

Therefore, we can logically expect the MSPs in the disc to be dif- 
ferent from the MSPs in the clusters. This is indeed the case. There 
is a greater proportion of single MSPs in the clusters and the ma- 
jority of the binary MSPs have very short orbital periods compared 
to those in the Galactic disc. In fact, many of the cluster binaries 
have properties similar to those of the rare eclipsing black widow 
pulsars seen in the Galactic disk population (King, Davies, & Beer 
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2003; Freire 2005). Recently, Bagchi & Ray (2009a, 2009b) have 
investigated the effect of different types of stellar collisions on the 
orbital parameters of the binary MSPs in globular clusters. They 
too find the cluster MSPs to be significantly different from those in 
the disc. 

In this work, we look at one of the important intrinsic parameters of 
the MSPs, namely the surface magnetic field. We find that a subset 
of slower MSPs in the globular clusters (for which field measure- 
ments are available) appear to have surface magnetic fields that are 
2-5 times larger compared to their disc counterparts. Even though 
the spin-periods of this subset of cluster pulsars are similar to the 
spin-periods of the MSPs in the disc. In order to understand this 
fact, we look at the spin-period and the magnetic field distributions 
of the MSPs in the disc and in the globular clusters and compare 
them to check for implied similarities/differences in their evolu- 
tionary histories. We also look at the details of the field evolution 
itself to see if the physics of the MSP formation in the disc and in 
the clusters could actually give rise to different surface fields. We 
examine whether close stellar encounters resulting in binary dis- 
ruption and/or re-formation of new binaries could be responsible 
for the higher magnetic fields observed in cluster MSPs. Or if this 
is an apparent effect caused by other external factors. 

To this end, first we discuss the MSP statistics in Sec. 2. In Sec. 3 we 
look at the details of the field evolution and the effect of the cluster 
conditions on the evolution. Finally, we summarize our conclusions 
in Sec.4. 



2 THE MILLISECOND PULSARS 

2.1 the definition 

It is well known that the MSPs occupy a distinctly separate region 
from the normal pulsars in the B — P plane. However, defining 
the MSP population with some accuracy is not simple since this 
requires separating out a class of neutron stars that have undergone 
a particular kind of binary evolution. The primary observed quanti- 
ties of a pulsar are P and P, and the most important derived quan- 
tity using both P and P is the dipolar component of the surface 
magnetic field given by (Manchester & Taylor 1977), 
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Figure 2. Histogram showing the distribution of spin-periods (P) of all 
known pulsars. The data correspond to the ATNF on-line catalog and Paulo 
Freire's on-line catalog of globular cluster pulsars, April 2010. 



possible for an MSP to have gone through very complex binary 
evolution. Therefore, in general, it is not easy to separate the MSPs 
from the normal pulsars using a definitive relation involving both 
P and P. 

In Fig. 1 we have plotted all known pulsars (with a reliable estimate 
for B) in the B — P plane, marking some of the B and P lines rel- 
evant for separating out the MSP population. Unfortunately, these 
simple criteria run into trouble with objects like the 16 ms pulsar 
J0537-6910 in LMC (Marshall et al. 1998). As this pulsar has a 
strong strong magnetic field {B ~ 10^^) and is also rather young 
(r ~ lO'^yrs), it is not likely to be a member of the generic recy- 
cled MSP population (this extra-Galactic pulsar is not included in 
Fig. 1). Therefore as a working hypothesis we shall use the classi- 
cal definition of P < 30 ms ( log(P/s) < — 1.5 ), barring such 
obvious misfits like J0537-6910. This way of defining MSPs does 
appear to have a certain merit. In Fig. 2 a histogram of the spin- 
period of all known pulsars is made. It appears that the pulsars with 
P < 30 ms may indeed represent a separate class as the period 
histogram shows a sharp dip around this value of P. 



Bs 



3.2 X 10^^(PP)^''^G 
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where P and P are in units of second and ss~ . Even though the 
term millisecond pulsar has traditionally been reserved for recycled 
pulsars with ultra- fast rotation (P< 30 ms) and a low magnetic field 
(B ~ 10* ~ 10^ G), the definition has mostly used the condition on 
P. Using P to classify pulsars is somewhat problematic because for 
a large number of pulsars (particularly for those in the clusters) a 
measurement of P is either not available or the value is not reliable 
due to contamination from the proper motion of the pulsar. 

The most reasonable criterion to define an MSP should be to use 
both P and P since the evolution of these two parameters are inter- 
related in the formation process of MSPs. Some effort has indeed 
been made to use such a relation assuming the MSPs to be recy- 
cled through straightforward LMXB evolution (Story, Gonthier, & 
Harding 2007). But the situation is different in globular clusters 
where, due to a much larger probability of stellar collisions, it is 



2.2 the statistics 

Interesting facts emerge when the MSP population, as defined 
above, is subjected to statistical analysis. To begin with, we sep- 
arate them into two groups - MSPs residing in the globular clusters 
(GC) and those in the Galactic disc (GD). Each of these groups 
are further divided into isolated and binary pulsars. Table. 1 shows 
the average values of P and B for each of these sub-populations. 
It needs to be noted that the number of pulsars with a measured 
value of B can be much smaller than the actual number of pulsars 
in a given sub-group (see the corresponding numbers in Table. 1). 
This is particularly true of the cluster pulsars. Because of this we 
have calculated another average of P using only those pulsars with 
a known B. This average, P^, is expected to be more commensu- 
rate with Bav than the Pav obtained using all the pulsars in the sub- 
population. The significance of these quantities shown in Table- 1 
are discussed below. 
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Figure 1. All known Galactic pulsars, with a reliable B measurement, in the B — P plane. The isolated and binary pulsars in the Galactic disc are marked by 
filled squares and open stars respectively. The isolated and binary pulsars in the globular clusters are marked by open circles and circles with dots within. The 
data is from the ATNF on-line catalog, taken in April 2010. 
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Figure 3. Kolmogorov-Smimov test for the P-distributions of different subgroups of the MSPs. The cumulative fractional distributions are shown for easy 
comparison. The solid, and the dashed lines correspond to the isolated and the binary MSPs in the globular clusters; whereas the dash-dotted and the dotted 
lines correspond to the isolated and the binary MSPs in the Galactic disc. The number of objects in each subgroup are shown within the parentheses. The K-S 
probabilities are indicated for each test. 



2.2.1 Spin period 

The well-known fact that the MSPs in globular clusters are, in 
general, spinning faster than their disc counterparts is immediately 
seen from the values in Table- 1. However, it can also be seen that 
the Pav for almost all the subgroups are similar except for the bi- 
nary MSPs in the disc which has a larger average spin-period than 
the rest of the subgroups. In other words, the disc binary MSPs are 



spinning more slowly than the rest of the MSPs. However, an av- 
erage difference of ~3 ms may not indicate anything significant as 
the standard deviation of P in these sub-groups range from ~ 4 - 
7 ms. Moreover, it is not clear if this could be due to any kind of 
observational bias. Though it should be noted that the size of the 
sub-groups are more or less similar except for the isolated objects 
in the disc (about a factor of 3 smaller). Therefore a bias (of any 
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Galactic Disc : 








isolated 


5.82(19) 


3.11 (18) 


5.97 (18) 


binary 


8.03 (54) 


9.32 (52) 


8.45 (52) 


all 


7.75 (73) 


7.72 (70) 


7.82 (70) 


Globular Clusters : 








isolated 


5.82 (61) 


14.70(11) 


6.16(11) 


binary 


5.58 (68) 


18.84(14) 


8.62 (14) 


all 


5.70(129) 


17.02 (25) 


7.54 (25) 



Table 1. The MSP Statistics. In the first column the type of the MSP 
sub-population has been indicated. The second and the third column show 
the average values of the measured spin-periods and derived dipolar field 
strengths respectively. P^^ in the fourth column is the average period cal- 
culated using only those pulsars for which B measurements exist. The spin- 
periods are in ms and the magnetic fields are in 10* G. The numbers in the 
parentheses denote the numbers of objects available for calculating the av- 
erages. The data is from the ATNF on-line pulsar catalog and Paulo Freire's 
catalog of globular cluster pulsars, April, 2010. 



kind) would be stronger for the isolated disc MSPs more tlian tfie 
rest of tfie groups. 

In order to understand tfie relative nature of the P-distributions 
of the various sub-populations, beyond the simple averages dis- 
cussed above we have performed the Kolmogorov-Smirnov tests 
on these populations. In Fig. 3 we compare the cumulative fractional 
P-distribution corresponding to each K-S test. And our conclusions 
from these tests are summarized below. 

(i) The isolated and the binary MSPs in the globular cluster itself 
are the least correlated {Pks ~ 0.002) and is in good agreement 
with that found by Hessels (2009). The result itself is surprising be- 
cause these two subgroups have very similar average spin-periods. 
Also the luminosity distributions of isolated and binary GC pulsars 
have been found to be statistically similar (Hessels et al. 2007). The 
low correlation is also contrary to expectations because all the pul- 
sars in the clusters are lilcely to have similar evolutionary histories, 
strongly influenced by stellar encounters. Since it is quite possi- 
ble for an isolated MSP to acquire a companion or a binary MSP 
to lose one given the high rates of stellar collisions in a globular 
cluster which effectively means that in a globular cluster the phase 
(isolated or binary) in which an MSP is observed at a given time 
could be quite temporary. 

(ii) The isolated and the binary MSPs in the galactic disc are 
also not correlated {Pks ~ 0.2) indicating major differences in 
their evolutionary histories. This is supported, to some extent, by 
the fact that their average spin periods are different (but not by any 
significant amount, as mentioned earlier). 

(iii) Similarly, there is little correlation {Pks ~ 0.03) between 
the spin-distributions of the binary populations in the disc and in the 
clusters. The explanation for this absence of correlation is possibly 
due to the fact that in the disc the MSP recycling typically happens 
in primordial binaries whereas in the clusters the binaries could 
well have formed by recent stellar encounters. Hence, the nature of 
the associated MSPs could be very different. 

(iv) Another surprising result is the very strong correlation 
{Pks ~ 0.9) between the the spin-distributions of isolated MSPs 
in the disc and in the clusters. The average spin-periods are simi- 



lar too, which may indicate similar evolution of the isolated MSPs 
in the disc and the clusters. But these two populations are not ex- 
pected to be so well correlated. In the disc the the isolated MSPs are 
supposed to evolve in primordial LMXBs and then evaporate their 
companions. In the clusters, this is certainly a possibility. However, 
the isolated MSPs are more likely to be results of stellar encoun- 
ters and consequent disruption of binaries in the clusters. It should 
also be remembered that the size of one of the subgroups (the disc 
pulsars) in this test is much smaller than the other and this could 
introduce some bias in the K-S test. 



2.2.2 Surface magnetic field 

First of all, it should be noted that the B measurement is available 
only for a small number of pulsars in the globular clusters (25 out of 
a total of 129). This is because of the difficulty in measuring the P 
there due to proper motion contamination. In fact, in many cases P 
is seen to be negative even when the system is non-accreting (and 
therefore has no reason to spin-up) indicating that the measured 
value is not at all reliable. In such cases there is no way to esti- 
mate B. Evidently, it would be relatively easier to measure larger 
Ps. This then automatically introduces a bias, that of preferentially 
measuring B in systems that have higher field values. On the other 
hand almost all the MSPs in the disc have a reliable B measure- 
ment. So any comparison between the disc and the cluster MSPs 
suffer from an inherent bias in this respect. 

Given the above problem it is still interesting to see that the mag- 
netic fields of the MSPs in clusters are larger by a factor of 2-5 
compared to their disc counterparts. Unlike in the case of P this 
difference is much larger than the standard deviation in B. Yet the 
corresponding Pj^s are similar in the disc and in the clusters, both 
for the isolated and the binary MSPs. However it should be noted 
that the fraction of MSPs, with field measurements, is much smaller 
in the clusters compared to the disc (see Tablel). For example, the 
fractions of MSPs with field measurements are 1 1/61 (isolated) and 
14/68 (binary) in the clusters. In contrast, the corresponding num- 
bers in the disc are 18/19 and 52/54 respectively. Consequently, in 
the disc Pav and P^^ are very similar for both the isolated and the 
binary pulsars, which, however, is not the case in the clusters. It 
can be seen that P^ is somewhat larger than Pav for the isolated 
cluster pulsars whereas P^JJ. is much larger than the corresponding 
Pav for the binaries indicating that B-measurements are probably 
being selectively made for slower pulsars in the clusters. Keeping 
this in mind, we can then conclude that the subset of slow cluster 
pulsars, for which B-measurements are available, appear to have 
higher magnetic fields compared to the disc pulsars that have on the 
average similar spin-periods to this subset. Therefore future mea- 
surements of the magnetic fields of the faster pulsars in the clusters 
are required to make proper comparison of these two sets. 

In Fig.4 we have plotted all the radio MSPs (as defined in Sec. 2) in 
the B-P plane, along with the AMXPs for which B measurements 
(mostly order of magnitude of upper limits) are available. It is seen 
that mostly, for a given value of P, the MSP with the highest value 
of B is a cluster pulsar. During the process of recycling a pulsar is 
spun up in the accretion phase. The maximum spin-up, for a given 
rate of accretion and a given strength of the surface field, is given by 
the following relation of spin-equilibrium (Alpar et al. 1982; Chen 
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Figure 4. All radio MSP candidates (left of the P = 30 ms line) in the B — P plane, along with the AMXPs (only estimated upper limits for B). The 
isolated and binary pulsars in the Galactic disc are marked by open and filled stars respectively. Whereas the isolated and binary pulsars in the globular clusters 
are marked by open circles and circles with dots within. The AMXPs are marked by open crosses. The line marked Spin-Up corresponds to Eq.2, defining 
the minimum spin-frequency attainable by recycling. The radio MSP data is the same as that in Fig. 1 . The AMXP data is taken from the references cited in 
Table-2. 
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Here M & R denote the mass and the radius of the neutron star, M 
& AfEd stand for the actual and the Eddington rate of mass accre- 
tion and B is the surface magnetic field. When M equals A^Ed the 
above relation defines the spin-up line i.e, the minimum period to 
which a pulsar with a given magnetic field can be spun-up. How- 
ever, even though all disc pulsars are consistent with the spin-up 
condition, some of the cluster pulsars like B 1821-24 or B 1820-30A 
appear to be above the spin-up line in the B — P plane. Given the 
errors and uncertainties in the B measurement, even if these pul- 
sars do not actually violate the above mentioned condition it can 
be said with some confidence that their field values are rather high, 
particularly in comparison with the disc MSPs with similar P. 

However, there is a serious caveat to this. Even among the disc 
pulsars, only a few (8 out of 18 isolated and 30 out of 52 binary 
MSPs with a B measurement) have P and hence B corrected for 
proper motions. And in those cases where such corrected values are 
available it is seen that the un-corrected values are systematically 
larger. For example, the average of the corrected B values for the 8 
isolated disc MSPs is 2.4 x 10* G whereas the uncorrected average 
is 2.8 X 10* G. The corresponding numbers in the case of binary 
disc MSPs are 4.7 x 10* G and 5.4 x 10* G respectively. Even 
though the errors, introduced in the P measurement by the proper 
motion contamination, are expected to be random it appears that in 
the case of the disc pulsars the B values uncorrected for the proper 
motion tend to be systematically larger than the corrected values. 
Admittedly, the difference is not very large for the disc pulsars. But 
in case of the clusters, where the proper motion is quite large, this 



difference could also be much larger. Therefore, this may make the 
measured field values in the clusters systematically larger. Still, it 
is not clear whether this effect is enough to explain the observed 
difference between the average field values of the disc and cluster 
pulsars. 

Once again, to understand the nature of the B-distributions better 
we have performed the K-S tests amongst the various subgroups of 
pulsars. However, except for the binary MSPs in the disc the sizes 
of the other subgroups are really small and hence the reliability of 
K-S test results are somewhat doubtful. We summarize our findings 
below and the corresponding plots have been shown in Fig.5. 

(i) The binary MSPs in the disc and the clusters have a low cor- 
relation {Pks ~ 0.2), as do the isolated and the binary MSPs in the 
disc (Pks^ 0.1). These results are in conformity of what we have 
seen for the respective P-distributions. And we can conclude that 
these subgroups have different evolutionary histories, as discussed 
before. 

(ii) Surprisingly, the B-distributions of the isolated pulsars in 
the disc and the clusters have extremely low correlation (Pks^ 5 x 
10"'*). This is in complete contrast with the high correlation 
(Pks ^ 0.9) of their P-distributions. But it should be remembered 
that for the cluster pulsars only a small number (11 out of 61) of 
objects have a measured B and hence only those have been used 
for this test. However, as mentioned earlier the Pb^s. of these two 
groups are quite similar, indicating selection of pulsars with similar 
P-values. 

(iii) Another surprise is the very good correlation (Pks ~ 0.7) 
between the isolated and the binary pulsars in the clusters, again in 
complete contrast with the poor correlation (Pks ~ 0.002) of their 
P-distributions. Though once again we need to note that a very 
small sample of objects (11 out of 61 for the isolated and 14 out 
of 68 for the binaries) are used for the B-test compared to the P- 
test. Whether this indicates a similar P-evolution for the groups is 
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Figure 5. Kolmogorov-Smimov test for the B-distributions of different subgroups of the MSPs, along with the cumulative fractional distributions. The solid, 
and the dashed lines correspond to the isolated and the binary MSPs in the globular clusters; whereas the dash-dotted and the dotted lines correspond to 
the isolated and the binary MSPs in the Galactic disc. The number of objects in each subgroup are shown within the parentheses. The K-S probabilities are 
indicated for each test. 
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IGRJ0029 1+5934" 


1.67 


< 3 X IQS 


Aql X-1 (1908+005)'' 


1.82 


<10^ 


XTEJ1751-305'= 


2.30 


3 - 7 X 108 


SAXJ1808.4-3658'' 


2.49 


1 - 5 X 10^ 


XTEJ1814-338<= 


3.18 


8 X 10® 


XTEJ0929-314^' 


5.41 


< 3 X 10^ 


SWIFT J1756.9-2508S 


5.49 


0.4 - 9 X 10® 


KS 1731-260'^ 


1.91 


< 7 X 10® 


EXO 0748-676' 


1.81 


~ 1 - 2 X 10^ 


*MXB 1730-335^ 


3.27 


~ 4 X 10* 



Table 2. AMXP and burst oscillation sources that have some estimate of 
their magnetic fields. MXB 1730-335 (oscillation frequency yet to be mod- 
ified) in NGC 6440 is the only one in a globular cluster with a field estimate 
The field values are inferred upper limits in most cases. The data is taken 
from a) Galloway et al. (2005), Burderi et al. (2007), Torres et al. (2008) 
b) Di Salvo & Burderi (2003), Casella et al. (2008); c) Wijnands et al 
(2005), Papitto et al. (2008); d) Di Salvo & Burderi (2003), Hartman et al 
(2008), Cackett et al. (2009); e) Papitto et al. (2007); f) Galloway et al 
(2002), Wijnands et al. (2005); g) Krimm et al. (2007), Pafi^no, Altami- 
rano, & Messenger (2010); h) Di Salvo & Burderi (2003); i) Loeb (2003) 
and j) Masetti et al. (2000) respectively. 



difficult to say because somewhat contrarily the Pg ^s are different 
for these two groups. 

In summary then we can think of the following possibilities giving 
rise to the higher surface fields observed in the cluster MSPs. 



(i) The higher field observed could simply be due to certain ob- 
servational biases (preferentially selecting high-B systems, mea- 
suring a higher value of P and hence B due to proper motion con- 
tamination and so on..). 



(ii) It is also possible that the cluster MSPs have higher fields 
because we are looking at a younger population in the cluster com- 
pared to the disc. These young MSPs observed today would slow 
down with time and migrate towards the right of the B ~ P plane. 
Looking at Fig.4 we can see that if the cluster MSPs move towards 
the right of the plot by appropriate amount then they may become 
identical with the current disc population. There is some support to 
this possibility from the observed AMXPs. The AMXPs turn into 
radio MSPs as soon as accretion stops in such systems. In princi- 
ple, their physical properties should be similar to very young radio 
MSPs. In table2 we list the AMXPs and the burst sources that have 
some estimate of their magnetic fields. These have been plotted in 
Fig.4 and we see that in the B — P plane the AMXPs occupy the 
same region as the cluster MSPs. 



(iii) On the other hand, the cluster MSPs could really have 
higher fields compared to the disc MSPs. This is possible only if 
the different nature of the recycling process in the clusters affects 
the field evolution significantly. In the next section we use a simple 
model for the evolution of the magnetic field to see if the stellar dy- 
namics in globular clusters facilitate halting the field decay earlier 
than is expected in a primordial LMXB. 
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3 RECYCLING : EVOLUTION OF THE MAGNETIC 
FIELD 
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In the standard formation scenario tlie MSPs are generated by 
recycling of ordinary pulsars in the LMXBs. Accretion-induced 
field decay is an integral part of this generic picture. A number 
of mechanisms, responsible for the evolution of the field, has been 
suggested. First, the magnetic field may be dissipated in the stel- 
lar crust by Ohmic decay, accelerated by heating as the accreted 
plasma impacts upon the star (Konar & Bhattacharya 1997; Urpin, 
Geppert, & Konenkov 1998; Brown & Bildsten 1998; Konar & 
Bhattacharya 1999a; Gumming, Arras, & Zweibel 2004). On the 
other hand, if the magnetic field resides in the superfluid core in 
the form of Abrikosov fluxoids then they may be dragged out of 
the core by the outward motion of superfluid vortices, as the star 
spins down (Srinivasan et al. 1990; Jahan Miri & Bhattacharya 
1994; Ruderman, Zhu, & Ghen 1998). This flux would then be 
subsequently dissipated in the accretion heated crust (Konar & 
Bhattacharya 1999b; Konenkov & Geppert 2001). Third, the mag- 
netic field could also be screened by accretion-induced currents 
within the crust (Bisnovatyi-Kogan & Komberg 1974; Lovelace, 
Romanova, & Bisnovatyi-Kogan 2005). In particular, the field may 
be buried under a mountain of accreted plasma channeled onto 
the magnetic poles (Hameury et al. 1983; Romani 1990; Brown 
& Bildsten 1998). When the accreted matter is large enough, the 
mountain spreads laterally, transporting the polar magnetic flux to- 
wards the equator and finally dissipating them there (Gumming, 
Zweibel, & Bildsten 2001; Melatos & Phinney 2001; Ghoudhuri 
& Konar 2002; Konar & Ghoudhuri 2004; Payne & Melatos 2004; 
Payne & Melatos 2007). 

In our earlier investigations we have assumed that the currents sup- 
porting the field finally get dissipated in the accretion heated crust, 
wherever they originally may have resided (Konar & Bhattacharya 
1997; Konar & Bhattacharya 1999a; Konar & Bhattacharya 1999b). 
We adopt the methodology developed in these articles (see Konar 
(1997) for details) for the present work and assume that purely 
crustal currents support the observed magnetic field of the neutron 
star. 

The evolution of the magnetic field is governed by the following 
equation 



riVK 

-=Vx(VxB) 



c^ 1 

— V X (- X V X B), 
47r o" 



(3) 



where bf V is the velocity of material movement and a is the elec- 
trical conductivity of the medium. The material movement in the 
crust due to mass transfer defines V. 

The spin-up of a neutron star, in a binary system, is caused by the 
angular momentum brought in by the accreted matter. In magne- 
tospheric accretion matter accretes with angular momentum spe- 
cific to the Alfven radius. Therefore, the total angular momentum 
brought in by accretion is : 



: SMRaVa, 



(4) 



where 5M is the total mass accreted, Ra and Va are the Alfven 
radius and Keplerian velocity at that radius. The final period of the 
neutron star then is : 



where Jns is the moment of inertia of the neutron star. The average 
dipolar field of the disc MSPs is ~ 10** G, whereas it is ~ lO'"' G 
for the MSPs in globular clusters (see Table- 1). So the minimum 
mass required to a spin a pulsar, with a surface magnetic field of 
~ 10^ G, to ~5 ms would be 5M ~ 0.01-0.1 Mq. The radius and 
the crustal mass of a neutron star remain effectively constant for 
this amount of accreted masses considered, and the corresponding 
change in the crustal density profile is negligible. We therefore take 
the mass flux to be the same throughout the crust, equal to its value 
at the surface. Assuming the mass flow to be spherically symmetric 
in the deep crustal layers of interest, one obtains the velocity of 
material movement to be 



V = 



M 



4:TTr^p(r) 



(6) 



where M is the rate of mass accretion and p(r) is the density as a 
function of radius r. 

The physical conditions of the crust enters Eq.3 through the elec- 
trical conductivity, a, which is a function of - i) the density of the 
current concentration, ii) the temperature of the crust and iii) the 
impurity content of the crustal lattice (negligible for an accretion 
heated crust). We have assumed the currents to be concentrated at 
a density ~ 10^"* gcm^'', and the impurity concentration to be 
~ 0.01 for a cold crust. These values are consistent with the fact 
that the magnetic fields of the isolated pulsars do not decay sig- 
nificantly over a time-scale of lO'' yrs (Bhattacharya et al. 1992; 
Hartman et al. 1997). 

The temperature of the crust is another important parameter for 
the field evolution. After the neutron star is born it monotonically 
cools down through copious emission of neutrinos (Page 1998). 
When, in the course of binary evolution, the neutron star actually 
starts accreting mass - the thermal behavior changes drastically. Ac- 
cretion releases heat and the crustal temperature quickly (in about 
^ 10^ yrs) settles down to a more or less uniform and steady value 
determined by the accretion rate (Miralda-Escude, Paczynski, & 
Haensel 1990; Zdunik et al. 1992). Though analytical expressions 
giving the crustal temperature for a given rate of accretion have 
been obtained by Zdunik et al. (1992) and Urpin & Geppert (1995), 
the values are too high for temperatures for M > 2 x 10^^" Mo/yr. 
Observations of thermal black-body radiation from the surface of 
accreting neutron stars indicate that for higher M the tempera- 
ture of the deep crustal layers probably saturate around ~ 10** K 
(see, for example. Brown & Gumming (2009)) . Accordingly for 
M > 10^^ Me/yr we assume the crustal temperatures to be in 
the range ~ 10* — 5 x 10* K. It should be noted that we consider 
only the temperature of the deep crustal layers, where the currents 
are assumed to be concentrated. And the temperature of the en- 
tire crust beyond a density of 10^" g cm^'^is practically the same 
though it drops by almost two orders of magnitude at the outer- 
most layers of the star (Gudmundsson, Pethick, & Epstein 1983; 
Potekhin, Ghabrier, & Yakovlev 1997). 

From the point of view of mass transfer (responsible for the phys- 
ical process of recycling), pulsars go through three distinct phases 
of evolution in typical LMXBs. These phases are as follows (Bhat- 
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tacharya & van den Heuvel 1991; Verbunt 1993; van den Heuvel 
1995) - 



ample the 'fly-by' kind of interaction between a single star and a 
binary). 



(i) The Isolated Phase - Though the stars are gravitationally 
bound, there is no mass transfer. In general, the isolated phase lasts 
between 10* — 10® years. 

(ii) The Wind Phase - The interaction is through the stellar wind 
of the companion which is likely to be in its main-sequence. In 
general this phase lasts for about 10*^ — lO** years with attendant 
rates of accretion ranging from about 10^^'' Mq yr^^to 10^^ M© 

(iii) The Roche-contact Phase - When the companion of the neu- 
tron star fills its Roche-lobe a phase of heavy mass transfer ensues. 
In this phase, the mass transfer rate could be as high as the Ed- 
dington rate (10^* M© yr^^for a 1.4 M© neutron star), lasting for 
< 10* years. But there has also been indications that the low-mass 
binaries may even spend ~ 10^" years in the Roche-contact phase 
with a sub-Eddington accretion rate (Hansen & Phinney 1998). For 
wide binaries, however, the contact phase may last as little as 10^ 
years. 



Earlier we have followed the evolution of the magnetic field for 
neutron stars in different types of binaries, assuming all the bina- 
ries to be primordial, undergoing standard phases of binary evolu- 
tion. For our adopted model (Konar & Bhattacharya 1997; Konar & 
Bhattacharya 1999a) it has been seen that the field decays rapidly 
in the initial phase, followed by a slow-down and a final freezing. 
The initial decay is due to the heating of the crust in which the 
currents undergo rapid ohmic dissipation. But as the accretion pro- 
ceeds there is addition of extra material in the outer layers of the 
star. Since an increase of mass makes a neutron star more com- 
pact, in the deeper layers of the crust this induces an inward radial 
motion. As a result the current carrying layers progressively move 
into higher density and higher conductivity region. Consequently 
the decay slows down. And when the entire current distribution, re- 
sponsible for the field, gets assimilated into the highly conducting 
(time-scale of diffusion larger than the Hubble time) core the decay 
stops altogether freezing the field at its final value. Understandably 
there is no further evolution of the field after freezing even if mass 
accretion continues. Also the higher the accretion rate the sooner 
the freezing sets in resulting in a higher value of the final surface 
field. 

In globular clusters, however, most binaries are not primordial. For 
example, the total observed number of LMXBs in globular clusters 
exceeds their formation rate in the disc by several orders of magni- 
tude, indicating a dynamical origin (Clark 1975). The composition 
of the binary itself may change, even more than once. Dynami- 
cal interactions in the clusters involving at least one neutron star 
are typically of the two and the three body types, as the probabil- 
ity of interactions involving four or more objects would be neg- 
ligibly small. Similarly, a three body interaction is essentially be- 
tween a single star and a binary system. Stellar interaction in the 
clusters involving a neutron star have been studied by a number 
of authors (Krolik, Meiksin, & Joss 1984; Rasio Sc Shapiro 1991; 
Davies, Benz, & Hills 1992; Davies & Hansen 1998; Rasio, Pfahl, 
& Rappaport 2000; King, Davies, & Beer 2003). We briefly discuss 
below the interactions which play an active role in pulsar recycling 
(see Camilo & Rasio 2005 for a detailed review). It needs to be 
mentioned that we do not consider the interactions that have no 
direct effect on the mass transfer phases of a binary pulsar (for ex- 



3.1 Two-Body Interactions 

The interaction could be a close tidal encounter or a direct phys- 
ical collision. However, the formation of stable binaries through 
tidal encounters is not very likely. Therefore, this process is not im- 
portant for MSP formation. For coUisional encounters there can be 
mainly two types of partners. 

Main- sequence star: Typically such a collision leads to the com- 
plete destruction of the main-sequence star forming a thick, rapidly 
rotating envelope around the pulsar. Depending on the amount of 
accretion the pulsar may be spun-up to millisecond periods or it 
may be only mildly recycled, either way giving rise to a single re- 
cycled pulsar. 

Red-giant star: Such a collision leads to the formation of a high- 
eccentricity binary. These collisions provide a natural formation 
process for eccentric low-mass binary pulsars with white-dwarf 
companions. If the post-collision neutron star - white dwarf bi- 
naries retain high eccentricities, then they could decay through 
gravitational-wave emission and possibly become ultra-compact X- 
ray binaries (UCXB) with Pb< 1 hr. These are important for pul- 
sar recycling since a number of the known AMXPs are actually 
members of UCXBs, which are also probably the progenitors of 
the black-widow MSPs (King, Davies, & Beer 2003). 



3.2 Three-Body Interactions 

When a single star interacts with a binary (the neutron star could 
either be a member of the binary or the single object) the result 
could be one of the following types of changes to the binary : 

(i) Exchange : In this kind of encounter one of the binary com- 
ponents is replaced by the single star. So a single pulsar could 
acquire a binary companion through this process. Alternatively, a 
previously formed binary pulsar could interact with another star or 
binary. This would lead to a new companion for an MSP, or for a 
non-recycled pulsar, or could release an MSP from a binary, cre- 
ating a single MSP. Systems with higher-mass companions, fast 
MSPs, and very high eccentricities, are likely to be the result of 
such exchange interactions, i.e., the presently observed companion 
was likely acquired later and is not the donor from which the pulsar 
was recycled. Exchange interactions between ordinary pulsars and 
primordial binaries also provide a natural way of forming possible 
progenitors of UCXBs (King, Davies, & Beer 2003). 

(ii) Disruption : Finally, it is possible for the binary to be com- 
pletely disrupted by its interaction with the single star. This would 
release a single pulsar. Depending on when the binary is disrupted 
the pulsar could be an MSP or a mildly recycled pulsar. 

(iii) Multiple Interactions : In some of the clusters the stellar 
densities are so high that the interaction time-scale for either ex- 
change or fly-by could be small making it possible for the binaries 
to undergo multiple interactions. Under the circumstances, a neu- 
tron star may actually go through various phases of evolution in a 
number of binaries sequentially. 
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Figure 6. Evolution of the surface magnetic field of a neutron star as a 
result of mass accretion, assuming the crustal currents to be concentrated 
at a density of 10^'* gcm^'^. Curves 1 to 5 correspond to mass accretion 
rates of M = lO"^^, 10"", 10"^", 10"'', 10"^ Mo/yr respectively. 
The calculations have been performed till the field decays no further. 



Given the above possibilities, a number of different situations can 
be envisaged where the phases of accretion are rather different from 
those in typical LMXBs. It has been seen that at least ~ 10"^ Mois 
required to spin a pulsar up-to millisecond periods. We have seen 
earlier that accretion of ~ 10"'^ Mois also enough for the field 
to freeze to its final stable value (Konar & Bhattacharya 1997). 
Because the mass of the crust of a typical 1.4M0 neutron star is 
~ 10"'^ Mq. And when this amount of material is accreted, the en- 
tire mass of the original crust containing the current carrying layers 
get accumulated into the highly conducting core where there can be 
no dissipation of the field. Once this amount of mass is accreted the 
field attains its final value even though the star may continue to go 
through further accretion and spin-up. Below, we list possible situ- 
ations where M> 10"^ Mq could be accreted onto a neutron star. 
For our calculations we have evolved the field till the final frozen 
value is attained. 

(i) The entire mass is accreted with M characteristic of the 
wind phase (M ~ 10"^'' — 10"^'^) in low-mass stars. This 
situation can arise if the original low-mass binary is disrupted 
before Roche-contact is established. The only way a pulsar could 
be spun up-to millisecond period entirely by wind accretion is if 
M ~ lO"^'^ — 10"^^ Mo/yr, assuming the main-sequence phase 
of an extremely low-mass companion to last for ~ 10^ — 10^" yrs. 



should be noted that collisions with red-giant stars are expected to 
give rise to UCXBs, extremely compact systems that are likely to 
have mass transfer at high M s. 

(iii) A brief phase of heavy mass transfer, followed by a long 
interval with low M , most of the mass being accreted in the 
second phase. This situation can arise either through an exchange 
interaction (the pulsar being always in the binary) or due to a 
binary (containing a pulsar) disruption followed by an acquisition 
of a new partner by the pulsar. We find that this is very similar to 
case A with the initial phase of heavy accretion having hardly any 
significant effect. 

(iv) A phase of heavy mass transfer followed by a brief interval 
of accretion with small M where most of the mass is accreted in 
the first phase. This is again very similar to case B. 



In summary, it can be said that given the time-scales of binary evo- 
lution and stellar collisions the amount of mass required to freeze 
the magnetic field to its final stable value would be mostly achieved 
in one particular phase of accretion with a given M. In Fig.6 we 
plot the evolution of the surface field with time for different values 
of M. It is seen that a higher magnetic field is retained by the pulsar 
if most of the mass is transferred in a short period with a high M 
even though initially the field decays faster due to a higher crustal 
temperature. This happens because the current carrying layers get 
assimilated into the core (which has effectively infinite conductiv- 
ity) faster for a higher M. A high rate of mass transfer can be re- 
alized if the pulsar is in a tight binary. Interestingly, it is already 
known that the binary MSPs in the clusters have relatively shorter 
orbital periods (tighter orbits) compared to the disc population. Dy- 
namically too, collisions of single pulsars with red-giant stars or ex- 
change interactions are likely to produce UCXBs - again systems 
where higher values of M can be achieved. So there definitely exist 
MSP formation channels in the clusters that are conducive of pro- 
ducing high magnetic field MSPs. But, in general, the average field 
values would be higher than the disc MSPs only if this were the 
dominant channel of MSP formation in the clusters. In an ongoing 
work we are looking at this question of relative importance of dif- 
ferent channels of MSP formation in the clusters (Bagchi & Konar, 
in preparation) but the results are not available yet. 



4 CONCLUSIONS 



(ii) A brief wind phase (or complete absence it), followed by 
heavy mass transfer (M ~ 10"^*^ — 10"* M©) characteristic of 
the Roche-contact phase in which most of actual mass accretion 
takes place. It needs to be mentioned here that this situation is very 
similar to typical LMXB evolution, except for a long wind phase 
with attendant low values of M prior to Roche-contact which is 
realized in most LMXBs. In our earlier investigations we have 
seen that the final field values attained with or without a phase of 
wind accretion are not greatly different (Konar & Bhattacharya 
1999a). Hence, even though the nature of binary evolution would 
be different the final P, B seen in the resultant MSP would be 
similar to one processed in a primordial low-mass binary. Now, 
this kind of mass transfer scenario is possible if an un-recycled 
(mildly or otherwise) pulsar acquires a partner which is already in 
an evolved phase (for example a red-giant star), in a tight binary. It 



In this work, we have compared the MSPs in the Galactic disc with 
those in the globular clusters vis- vis their distribution of the spin- 
period and the surface magnetic field. The statistical nature of the 
two populations are as follows. 

(i) The average spin-periods of isolated and binary MSPS in 
globular clusters as well as isolated MSPs in the Galactic disc are 
very similar. Though the average spin-period of the binary MSPs in 
the disc is somewhat different, the difference is not very significant. 

(ii) The Kolmogorov-Smimov probabilities indicating the ex- 
tent of correlation of the spin-period between various subgroups 
of MSPs are found to be as follows - 

• isolated vs. binary in the clusters - Pks ~ 10""^, 

• isolated vs. binary in the disc - Pks ~ 10"^, 
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• isolated in the clusters vs. isolated in the disc - Pks ~ 0.9, 

• binary in the clusters vs. binary in the disc - Pks ~ lO^'^. 

(iii) Finally, the subset of cluster MSPs, for which field mea- 
surements are available, appear to have 2-5 times higher surface 
magnetic fields compared to the disc MSPs. Though the MSPs 
with field measurements are actuallay a slower subset of the clus- 
ter MSPs their average spin periods are similar to the MSPs in the 
disc which are, in general, slower than the cluster MSPs. We feel 
the difference in the surface magnetic field could be due to one or 
several of the possibilities listed below. 

• There are systematic biases in P (hence B) measurement 
for cluster MSPs. 

• The cluster MSPs are younger and may evolve to a distri- 
bution similar to the disc MSPs with time. 

• Preferential recycling of MSPs in tighter binaries with high 
rates of attendant mass transfer may actually result in cluster pul- 
sars retaining higher magnetic fields. 
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